Prepared for submission to JCAP 



Gravitational Waves in Ghost Free 
Bimetric Gravity 



Morteza Mohseni 

Physics Department, Payame Noor University 19395-3697 Tehran, Iran 
E-mail: m-mohseni@pnu.ac.ir 

7^ , ' Abstract. We obtain a set of exact gravitational wave solutions for the ghost free bimetric 

theory of gravity. With a flat reference metric, the theory admits the vacuum Brinkmann 
plane wave solution for suitable choices of the coefficients of different terms in the interaction 
potential. An exact gravitational wave solution corresponding to a massive scalar mode is also 
admitted for arbitrary choice of the coefficients with the reference metric being proportional 
to the spacetime metric. The proportionality factor and the speed of the wave are calculated 
in terms of the parameters of the theory. We also show that a F(R) extension of the theory 
admits similar solutions but in general is plagued with ghost instabilities. 
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1 Introduction 

Following the pioneering work by Fierz and Pauli [1], there have been numerous efforts to 
construct a consistent theory of massive gravitons, see [2] and [3] for reviews. Most of the 
resulted models suffered from subtle theoretical problems such as the van Dam-Veltman- 
Zakharov discontinuity [4, 5] and the appearance of ghosts and subsequent instability [6], 
which make them far from being conclusive. Such efforts have recently been intensified, partly 
as a consequence of the belief that massive gravitons might explain the observed accelerated 
expansion of the universe [7-18]. In a remarkable breakthrough, a nonlinear extension of the 
Fierz- Pauli model has been introduced recently [19, 20] which is shown to be free of the 
Boulware-Deser ghosts [21]. Thus, the model is usually referred to as the ghost free massive 
gravity. This model has attracted a lot of attentions, a nonexhaustive list of the relevant 
references includes [22-33]. 

A particular feature of the massive gravity theories is the simultaneous existence of at 
least two metrics, the dynamical or full metric which specify the geometry of the spacetime, 
and the absolute or reference metric which is related to the graviton mass term. In the ghost 
free model of massive gravity, the reference metric is fixed while the full metric is dynamical, 
in general. On the other hand, in Ref. [34], it has been shown that considering a dynamical 
reference metric does not destroy the ghost free property of this nonlinear theory of massive 
gravity. Thus, it would be interesting to construct a model in which both metrics are dynam- 
ical. Actually, this has been done in Ref. [35] in the context of a bimetric theory. Basically, 
this is an extension of the ghost free model of massive gravity in which the reference metric 
has its own kinetic term (an Einstein-Hilbert-like term) and an interaction term between the 
two metrics with a form similar to the nonlinear potential of the ghost free model of massive 
gravity. 1 Spherically symmetric solutions of the model have been studied in Ref. [37]. The 
cosmological solutions to the model have been presented in Refs. [38] and [39]. A mechanism 
which removes the ghost from the model is discussed in Ref. [40]. An equivalent model has 
been formulated in terms of vielbeins in Ref. [41]. A modified version of the model has been 

1 Actually, it seems that there is no standard terminology regarding these two types of metrics. In Ref. 
[36], the terms "background" and "foreground" are used for these metrics. Here, we reserve the term reference 
for the metric introduced to define the graviton mass and the term spacetime for the one which determines 
the causal structure of spacetime. 
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constructed [42] in which the curvature scalar of the spacetime metric is replaced with an 
arbitrary function of it using an the auxiliary scalar field and it was shown that the modified 
model admits a cosmological solution corresponding to an accelerating universe. In Ref. [43] 
it was shown that the model can be extended by replacing the curvature scalars with the 
Lovelock invariants. 

In this work we seek classes of exact pp-wave solutions of the ghost free bimetric gravity. 
The first motivation for this comes from the general interest on gravitational waves. In 
fact, gravitational waves are among the most important theoretical predictions of theories 
of gravity and are expected to be observed in near future. Thus, one expects any modified 
theory of gravity to accommodate such solutions. They are also of great interest in other 
areas of theoretical physics such as string theory, see, e.g., [44] and references therein. Here, 
we present a set of massless and massive exact gravitational pp-wave solutions. This is 
particularly interesting because most of the known gravitational wave solutions to modified 
theories of gravity, models of massive gravity in particular, are obtained within the framework 
of linearized approximation, see e.g., [45] and references therein. For the case of massive 
gravity, which is not free of ghosts at the linear order, it would be more plausible to search 
for solutions describing strong gravitational waves as well. An Aichelburg-Sexl-type plane 
wave solution has been found for the ghost free massive gravity [27] and a few other modified 
theories of gravity [46] . The second motivation lies in the bigravity model itself. Investigating 
new solutions to the model sheds more light on the physical content of the model and its 
potential strengths. The data associated with new solutions, including gravitational wave 
ones, might be used to constraint the model and its parameters. 

In the following sections we start with a review of the ghost free model of bimetric gravity 
and then present massless and massive pp-wave solutions in the absence of matter fields. We 
show that the model admits vacuum plane gravitational solution for a suitable choice of the 
parameters by choosing a flat reference metric. For massive solutions to be admitted, we can 
choose a reference metric proportional to the spacetime metric. We show that the speed of 
the wave can be calculated in terms of the parameters of the theory which at the same time 
fix the proportionality factor. We also show that the F(R) extension of the theory admits 
these solutions but it suffers from ghost instabilities. 

2 The ghost free bigravity 

The action for the ghost free bimetric gravity [35] may be written in the following form [39] 



where R and 7Z are the curvature scalars associated with the spacetime metric g^ u and the 
reference metric fn V respectively, C m is the Lagrangian density for the matter fields generically 
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is the interaction potential for the two metrics g^ u and f^. Here, /3 n are arbitrary constants 
and 



eo = 1, 
ei = [X], 

e 2 = \ {[X? - [X 2 ]) , 

e 3 = i([X] 3 -3[X][X 2 ] + 2pr 3 ]), 



(2.3) 



e4 = det(X) 

where square brackets stand for the trace with respect to the spacetime metric g^ u . Thus, 
the parameters of the model are the five constants /3 n , the coupling constants M p and M p , 
and the mass m. Variation of the above action with respect to g^ u and results in the 
following equations of motions, respectively 



2 3 i 



n=0 

f 2 



m 2 M 2 3 



7 /XI/ 
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V n=0 



£(-i) n /9 4 _„ f/ M Y ( ^(x') + LxY^(x'\ 



(2.4) 
(2.5) 



where is the Einstein tensor for g^ Ul Qn V is the Einstein tensor associated with f^ u , T^ u 
represents the energy-momentum tensor, X' = \J / and Y/ n \(X) are given by 



Y {Q) (X) = 1, 
Y (1) (X) = X - [X], 

Y (2) (X) = X 2 -X[X] + ^([X] 2 -[X 2 }), 

Y (3) (X) = X 3 - X 2 [X] + \x ([X] 2 - [X 2 ]) - I ([Xf - 3[X][X 2 } + 2[X 3 ]) . 



(2.6) 



3 Massless gravitational waves 

In this section we look for massless plane gravitational wave solutions propagating in a flat 
background. Thus, adopting the chart (u, v, x, y) in which u = t — z and v = t + z are the 
null coordinates, we set 

/ -\ 0\ 



which is the Minkowski metric, and 
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where H = H(u, x, y) is the wave profile. This corresponds to the following line element 



ds 2 = —dudv — H(u, x, y)du 2 + dx 2 + dy 2 
which is the Brinkmann expression for the plane wave. This results in 



(3.3) 
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For the positive roots, inserting these into Eq. (2.4) results in 
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from which we obtain 
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Similarly, the other equation of motion, Eq. (2.5), reduces to 
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into account, this results in 



in which a, 



r/3 n , and V 2 is the Laplacian in the transverse plane. Taking Eq. (3.5) 
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V 2 H 
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which upon substitution of the above obtained values of 02, 3, 4 reduces to 

V 2 H = 



(3.6) 
(3.7) 

(3.8) 



which is the standard equation for the wave profile in general relativity. For the negative 
roots in Eq. (3.4) we obtain the same results but with (3\ — > — f3\ and 03 —> —03- 

Here, two more points are in order. First, because in the absence of matter fields the 
action is invariant under the following interchanges [39] 



9iw f^u, Pn 04-n, M p O M p , m 2 



m 2 M 2 



(3.9) 
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we have automatically another solution in which f^ v is equal to the Minkowski metric and 
is the plane wave. Second, if we generalize the model by replacing 1Z in Eq. (2.1) with an 
arbitrary function of it, i.e., J-(1Z), then the first term in the left-hand side of Eq. (2.5) gets 
modified as 

dn 'V - + (V u (/) - y u)^(f)v) dn 

while the other terms (and also the other equation of motion, Eq. (2.4)) remains intact. 
It can easily be shown that it also admits the same solution as presented above provided 
J-(0) = 0. In the absence of matter fields, the same thing is also true for the case where the 
curvature scalar of the spacetime metric is promoted to F(R), as in Ref. [42]. It should be 
noted that in general, the extra degrees of freedom associated with these extensions may lead 
to appearance of ghosts in the model. In some specific situations, these ghosts can be avoided 
[42]. 

4 Massive gravitational waves 

In this section we seek massive pp-wave solutions in vacuum. We switch to the Cartesian 
coordinates (t, x, y, z) and start with the following ansatz 

= K(fit - z)r,^, (4.1) 
V = a 2 K{(3t - z)r,^ (4.2) 

in which < (3 < 1 and a are constants, and the specific form of the function K{f3t — z) is to 
be determined from the equations of motion. With this ansatz, we have 

X% = ±a5£ (4.3) 

and 

X% = ±a- 1 5» (4.4) 
from which we obtain (for the positive signs) 
3 

Y,(- l TPn{9^ n)v {X) + g vX Yf n)ll {X)) = 2(A) + 3afr + 3a 2 /3 2 

n=0 

+a 3 p 3 )K(pt - z) Vll „ (4.5) 

3 

^2(-l) n Pn(UxY^ u (X>) + UxYfa^X')) = 2(a 2 /3 4 + 3a/3 3 + 3/3 2 



71=0 



+a- 1 p 1 )K(l3t - z)r,^. (4.6) 



On the other hand, the only nonvanishing off-diagonal component of the Einstein tensor 
associated with the metric given in Eq. (4.1) is 

2K(/3t - z)K"(/3t -z)- 3(K'(pt - z)) 2 
* z "2 (K(pt-z))* [ } 

where primes denote differentiation with respect to fit — z. Now, Eq. (2.4) requires Gt z to 
vanish and hence 

K(/3t-z) = K {fit-z)~ 2 , (4.8) 
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with Kq being a constant. Now, the Einstein tensor reduces to 

G Mi , = -3(l-/3 2 )(/3t-z)-V, (4.9) 
and a similar relation for Q^. Inserting these into Eq. (2.4) we obtain 

2 

1 - /3 2 = ^-(A) + 3afr + 3a 2 /3 2 + a 3 p 3 )K . (4.10) 



n fit IVl^ H 

P = ~~o1Ta2~ ( a 04 + 3a/3 3 + 3/3 2 + a" 1 /?!)^- (4.11) 



Similarly, from Eq. (2.5) we obtain 

m 2 M? 
~3M p 

Thus, the model admits the solution given in Eqs. (4.1) and (4.2) provided the above two 
conditions holds. This corresponds to a scalar massive wave propagating in the z- direction 
with subluminal speed j3. The above two equations determine the values of /3 and a in terms 
of the parameters of the model. As an example, by choosing /3 2 to be the only nonvanishing 

/3 n , one obtains a = -j^- and /3 = Wl Then the constraint j3 < 1 results in 

M 2 

< K0P2 < m 2^ji ■ Starting with the negative roots in Eqs. (4.3) and (4.4) results in similar 
relations but with f3\ — > —f3\ and 03 —> —fa- 

5 The F(R) extension 

For the extended model where in action (2.1) the curvature scalar R is replaced with a generic 
function F(R), we obtain, instead of Eq. (2.4), the following equation of motion 

1 23 

(iV + g„ v U - y,V v )F\R) - -g, u F{R) = ^(-l) n /3„ (g^Y^X) + /x O v) , 

n=0 

(5.1) 

while the other equation of motion, Eq. (2.5) remains unchanged. Here, F'(R) = dF d ^ 
and we have assumed that matter fields are absent. Now, repeating the calculations of the 
previous section, one can show that the massive pp-wave solution given in Eqs. (4.1) and 
(4.2) is still admitted provided Eq. (4.10) is replaced with 

1 - /3 2 = o "f„ , (fa + 3a/3i + 3a 2 /3 2 + a 3 /3 3 + ^^-) K , (5.2) 



3F'(R ) V 2m 2 

where the curvature scalar R$ = 12(1 — (3 2 )Kq 1 is constant. The consistency of the equations 
constraints the function F(R). For the example case where /3o, 1,3,4 = 0, /3 2 7^ 0, we obtain 

a = ^F'(Rq) (J^ ~ 6m^F^ ]fe) which re q uires that the relation < fe holds. 

The massive wave solution is a constant curvature spacetime which can be used to 
explore the properties of the F{R) extension further. This can be achieved by considering the 
propagating modes resulted from the extension. Taking the trace in both sides of Eq. (5.1), 
we obtain 

3DF'{R) + RF'(R) - 2F(R) = A (5.3) 
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where A = - m ^ 3 n=0 (-l) n P n Y^(X). Now, defining = F'(R) and % = \{2F{R) - 
RF'(R) + A), we reach the following Klein-Gordon equation 

dV 

n *=* (5 - 4 » 

in which V plays the role of an effective potential. For the spacetime under consideration, 
the minimum of V occurs at <j)Q where the following relation is satisfied 

2F(R ) - RoF'(Ro) + A = 0. (5.5) 

For the perturbation eft = cfto + 6(f) around the minimum value, we have 

(□ - Ml R )5cp = (5.6) 

where M FR = | ( F"(Ro) ~ ^oV Since one expects F"(Rq) to be very small [47], Mp R 
is positive. This gives rise to a massive propagating mode essentially similar to the F{R) 
extension of massless gravity [48]. However, the F(R) extension also excites ghost modes in 
the model. To see this more clearly, we first rewrite the action in terms of g^ u = K(j3t — z)r]^ u 
and /„„ = L{f3t — z)r]^ u with the assumption that matter fields are absent. This reads 
S = J Sd 4 x where 



2 V 2 7 I K 2 V 2 7 V L J 

+m 2 M 2 ((3 K 2 + APxKVKL + GfoKL + 4f3 3 LVKL + f3 4 L 2 ^ (5.7) 

where primes stand for differentiation with respect to u = fit — z. The associated equation of 
motion of K(j3t — z) can be obtained from the Euler-Lagrange equation 

EM^fs^i-O (5.8) 



which gives 



du n \ d, 

n=0 \ V du n 



3, ^ ( 2K" K' 2 \ 9 f„ „ fh „ L „ L 




4 (1 - P Z ) ~ -^j = 2m 2 + 3^y- + 3A- + fc-\f 

and a similar equation for L(f3t — z). These equation have consistent solutions K(f3t — z) = 
K~o(f3t — z)~ 2 and L(f3t — z) = a 2 K(/3t — z) which upon insertion in Eq. (5.9) results in Eqs. 
(4.10) and (4.11). 

For the extended model in which the scalar curvature R (in the spacetime sector or 
the reference sector or both) is replaced with a generic function F(R), the relevant kinetic 
terms in the right-hand side of Eq. (5.7) should be modified accordingly. For the case of the 
spacetime sector, the first term on the right-hand side of Eq. (5.7) should be replaced with 



M 2 I 3 



2K" K' 2 



K 2 K 3 
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With this, the first two terms in the left-hand side of Eq. (5.8) are as follows 



d? ( dF(R) \ d f dF(R) K'\ d?F(R) „ d*F(R) 2 
du 2 V dR ) du\ dR K ) dR 2 dR 3 [ ' 

K' d 2 F(R) 
+ ~K dR 2 

M 2 

where we have dropped overall factors — j- and 3(1 — /3 2 ). In general, noting the relation 

R = 2 j^ 1 — > the first three terms in the right-hand side of the above relation contain K"" 
and K'". Thus, the equation of motion for K will be of the following general form 




d 4 K(u) n d 3 K(u) „d 2 K{ 



a 



Thus, by considering a perturbed solution 

K = + 5K 

where = K (/3t — z) 2 , we will end up with an equation of motion for 5K which is 

generally of fourth order derivatives. According to the Ostrogradski theorem, such equations 
result in ghost propagating modes [49]. Therefore, the F(R) extended model is plagued with 
ghost instabilities. Similar arguments hold when the extension is made in the reference sector. 



6 Conclusions 



We showed that the ghost free bimetric theory of gravity admits exact gravitational pp- 
wave solutions. For vacuum, the Brinkmann gravitational wave solution is admitted with 
a flat reference metric for suitable choices of the parameters. The resulting constraints on 
the coefficients of different terms in the potential determines three of them in terms of the 
remaining ones, namely, those of the zero and first order terms. We showed that a F(R) 
modified version of the theory admits the same solution. An exact massive gravitational 
wave solution is also admitted with a reference metric being proportional to the spacetime 
metric. The proportionality constant and the speed of the wave were calculated in terms of 
parameters of the model. The massive solution corresponds to a scalar mode of gravitational 
wave. We showed that the F(R) extension of the model also admits the solution but the 
massive solution is not stable because of the presence of ghost modes. The exact gravitational 
wave solutions presented here are useful for studying the strong gravitational waves in the 
context of ghost free bimetric gravities and are more compatible with the nonlinear nature of 
such theories. 
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